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In common descriptions of phase transitions, first order transitions are characterized by discon¬ 
tinuous jumps in the order parameter and normal fluctuations, while second order transitions are 
associated with no jumps and anomalous fluctuations. Outside this paradigm are systems exhibiting 
‘mixed order transitions’ displaying a mixture of these characteristics. When the jump is maximal 
and the fluctuations range over the entire range of allowed values, the behavior has been coined an 
‘extreme Thouless effect’. Here, we report findings of such a phenomenon, in the context of dy¬ 
namic, social networks. Defined by minimal rules of evolution, it describes a population of extreme 
introverts and extroverts, who prefer to have contacts with, respectively, no one or everyone. From 
the dynamics, we derive an exact distribution of microstates in the stationary state. With only two 
control parameters, Nps (the number of each subgroup), we study collective variables of interest, 
e.g., X, the total number of I-E links and the degree distributions. Using simulations and mean-held 
theory, we provide evidence that this system displays an extreme Thouless effect. Specihcally, the 
fraction X/ {NjNe) jumps from 0 to 1 (in the thermodynamic limit) when Ni crosses Ne, while all 
values appear with equal probability at Nj = Ne- 


I. INTRODUCTION 

In systems with many interacting degrees of free¬ 
dom, interesting collective phenomena are associated 
with phase transitions, e.g., in ferromagnetism. Here, a 
suitable macroscopic variable characterizing the state of 
the system - the order parameter - typically changes its 
behavior in some rather dramatic fashion. In standard 
textbooks, phase transitions are classified by the Ehren- 
fest scheme: first order, second order, etc. We also learn 
to expect certain characteristics associated with each or¬ 
der. Thus, across the first order transition, the order 
parameter jumps, while its fluctuations are ‘normal’ (on 
either side). Metastability, hysteresis, and co-existence 
are other common features associated with this kind of 
transition. By contrast, opposite characteristics, e.g., no 
discontinuity and anomalously large fluctuations, are as¬ 
sociated with second order transitions. 

Though such properties are observed in most physical 
systems, there are exceptions. In the context of one¬ 
dimensional Ising models with long range interactions, 
Thouless EHi found ‘mixed order’ transitions, at which 
the order parameter jumps discontinuously and exhibits 
large fluctuations. Since then, several systems with such 
properties have been discovered |4l4l0l|. In particular, 
the term ‘extreme Thouless effect’ was coined recently 
0[ll to describe a case where, at the transition, both 
the jump and the fluctuations are maximal. In this paper, 
we report another system displaying such an effect, in 
the context of a minimal model of social interactions, 
involving dynamic networks with preferred degrees. 

In our previous studies of such networks , we 

introduced the notion that an individual (i.e., node) 


adds/cuts links to others according to its ‘preferred de¬ 
gree’: K. The evolution of the simplest version of such 
networks is: In each time step, a random node is cho¬ 
sen and its degree, k, is noted, li k > k, the node cuts 
one of its existing links at random. Otherwise, it adds a 
link to a randomly chosen node not connected to it. In 
the steady state of a homogeneous system (all nodes as¬ 
signed the same k), this ensemble of apparently random 
graphs displays quite different properties [l5l| from the 
standard Erdds-Renyi case [l^. Taking a small step to¬ 
wards describing an inhomogeneous society, we consider a 
heterogeneous system of two subgroups, with Ni ^2 nodes 
assigned different /t’s. Letting ki < K 2 , we naturally refer 
to the first group as ‘introverts’ (I) and the latter one as 
‘extroverts’ (E). Despite the simplicity of its rules, such 
a system exhibits a rich variety of properties, discovered 
mainly through simulations [l4| . On the analytic front, 
progress has been modest, since the underlying dynam¬ 
ics violates detailed balance and the stationary state will 
have non-trivial persistent probability currents 0. Un¬ 
der these circumstances, to gain some insight, we turn 
to limiting cases which embody the main features of the 
full system. In this spirit, we consider the ultimate limit: 
Ki = 0, K 2 = oo. In other words, these are extreme 
introverts and extroverts (or XIE, for short). Specifi¬ 
cally, we are able to find an exact expression for the sta¬ 
tionary distribution, obtain analytic predictions for var¬ 
ious quantities which are confirmed by simulations, and 
provide good evidence that the transition across Ni = N 2 
displays an extreme Thouless effect. While preliminary 
results have been reported earlier Bin, we will present 
a more detailed study of this model here. 

First, for the readers’ convenience, we provide a sum- 
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mary of the preliminary findings, by including a complete 
description of the model, the master equation associated 
with the stochastic process, the exact microscopic dis¬ 
tribution in the steady state, and a mapping to an two- 
dimensional Ising model with peculiar interactions. In 
Section III, we study the statistical properties of X (the 
total number of I-E links), an ‘order parameter’ which 
corresponds to the magnetisation of the Ising model. Us¬ 
ing Monte Carlo techniques, we report findings much be¬ 
yond those in ref. |1§ . [l9l| : including a power spectrum 
study of the time traces of X, as well as some first steps 
towards a finite size analysis for X, as a function of Xj^e 
(the number of I,E’s). These provide further evidence 
for the principal characteristics associated with an ex¬ 
treme Thouless effect. The following section is devoted 
to investigations of a standard characterization of net¬ 
works: degree distributions. In the Ising language, these 
correspond to novel measures of the system, offering a 
more detailed picture than the magnetisation. For sys¬ 
tems with Nf yf Ne, predictions (i.e., no fitting param¬ 
eters) of a self-consistent mean-field theory are largely 
confirmed by simulations. We end with a summary and 
outlook, while the Appendices contain many of the tech¬ 
nical details. 


II. THE XIE MODEL AND THE STATIONARY 
DISTRIBUTION 

With the motivations for this model presented in both 
the introduction and previous studies [l^, [I^, we pro¬ 
vide the specifications of our model, using the language 
of a social network. Our population consists of N indi¬ 
viduals, divided into two groups: Nj introverts and Ne 
extroverts. Their behavior is ‘extreme,’ in the sense that 
the former/latter prefers contacts to none/everyone. The 
rules of evolution cannot be simpler: 

• In each time step, a random individual is chosen. 

• If an introvert is chosen, it cuts a random existing 
link. 

• If an extrovert is picked, it adds a link to a random 
individual not already connected to it. 

Note that one link changes at every step, except when 
the chosen individual is ‘content,’ i.e., a totally isolated 
/ or a fully connected E. Obviously, such a contented 
pair cannot be present simultaneously in our model. In 
this sense, our system may be termed ‘maximally 
frustrated,’ a measure |14| we will not pursue here. 

In simulation studies, it is customary to define one 
Monte Carlo Step (MCS) as N such attempts, so that 
each node has an even chance of being chosen after N 
attempts. Our main interest will be the statistical prop¬ 
erties of this system at long times, once it settles into 
the stationary state. We emphasize again that, in this 



FIG. 1. (Color online) The nodes of the two groups are de¬ 
noted by circles: blue open (!) and red closed (E). The black 
lines represent the active cross-links and the red dashed lines, 
the frozen E-E links. For this network, the sets of fc’s are: 
ki = {1, 0, 1,1,1, 2,1}, and Re = {6, 5, 4, 4}. Thus, this con¬ 
figuration contributes 1, 5,1 to pi {k = 0, 1, 2) and 2,1,1 to 
pE {k = 4, 5, 6), respectively. 

minimal model, there are just two control parameters: 
Nj^e- For large Y, we may consider different ‘thermo¬ 
dynamic’ limits, e.g., fixed difference Ne — Nj or ratio 
Ne/Nj. Note also that, though the total number of links 
can reach N {N — 1) /2, the maximum number of cross¬ 
links between the two groups is 

J\f = NiNe (1) 

Clearly, regardless of how the system is initialized, all 
the intra-communities links (/-/ or E-E) will quickly be¬ 
come static (all absent or present). Only the I-E cross¬ 
links are dynamic, depending on which node happens to 
be chosen (illustrated in Fig. [1}. In other words, we 
may limit our attention to the set of bipartite graphs. 
Of course, the total number of cross-links, X (g [0,A/]), 
fluctuates and should display interesting statistical prop¬ 
erties. The first impression of such a minimal model is 
that it must be trivial. In particular, it may be argued 
that, since the probability for a link to be cut or added 
is proportional to Nj/N or Ne/N, the fraction 

X = X/N (2) 

should be simply Ne/N. In the rest of this article, we 
will show the behavior of x to be dramatically different. 

The configuration space here consists of the Nj x Ne 
incidence matrices, N. We denote its elements by riy 
(f G [1, Af/] ,j G [1, A^)]) which is 1 or 0 when the link 
between an introvert node i and an extrovert j is present 
or absent, respectively. Clearly, this configuration space 
is identical to that of a Nj x Ne Ising model on a square 
lattice. Following the lattice gas language [l^, we will 
refer to n = 1,0 as a ‘particle’ or a ‘hole.’ Meanwhile, 
the dynamics of cutting/adding corresponds to a kinetic 
Ising model with spin flip dynamics [2ll| (i.e., without 
particle conservation). Since X = JiijUij and Af is the 
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total number of sites in the Ising system, we have the 
mapping 


a; = (1 + m) /2 (3) 

where m G [—1,1] is the magnetisation. To further this 
correspondence, let us define 

h = A/TV; A = Ne-Ni (4) 

with h playing the role of the magnetic field. Thus, our 
interest in how x responds to Nj^e translates into finding 
a ‘equation of state’ m{h). In this language, the naive 
expectation is trivial (and false): m (h) = h. 

Of course, unlike in the Ising case, the statistical me¬ 
chanics of the XIE model is defined by dynamical rules 
rather than a Hamiltonian. Therefore, to study quanti¬ 
ties of interest, we must first find the distribution of the 
stationary state, (N), as opposed to just writing a 
Boltzmann factor. 

With the dynamics specified, we can write the master 
equation for P(N, t), the probability of finding the system 
in configuration N at time t: 

V{N,t + l)-V{N,t) = 

[W(N, N')V{N', t) - W(N', N)V{N, t)] (5) 

{N'l 

where IT(N',N) is the probability for configuration N to 
become N' in a step (an attempt). For W, we note that 
the dynamics for an I node i involves cutting a random 
existing link, so that ki, the number of links it has (i.e., 
particles in row i in N) will be needed. Similarly, for a E 
node /, we will need pj, the number of holes in column j 
in N. Letting hy = 1 — Uij, we define 

ki = YjjTlij^ Pj = YjiTlij ( 6 ) 

Clearly, these variables reveal a not-so-explicit symme¬ 
try in the dynamics of XIE. Similar to the Ising spin 
flip symmetry (n n), there is an additional, transpose 
operation: 


Uij Uji 0 TV/ Ne (7) 


We will refer to this, a la Ising, as ‘particle-hole sym¬ 
metry,’ which will play an important role in discussions 
below. A layman’s way to phrase this symmetry is: The 
presence of a link is as intolerable to an introvert as the 
presence of a ‘hole’ is to an extrovert. 

With these preliminaries, the IF(N',N) in Eqn. ([S]) 
reads 


E 


e(fcO 

h 


n'.^n^j 


Q(Pj) 

Pj 


n'.jUij 


TV 


( 8 ) 

where 0 (x) is the Heavyside function (i.e., 1 if x > 0 and 
0 if X < 0) and the product of S’s ensures that only one 
Uij may change in a step. It is straightforward to verify 
that Eqn. ([S]) respects particle-hole symmetry. 


The dynamics defined by Eqn. m is clearly ergodic. 
More remarkably, unlike those in less extreme models 
of introverts and extroverts [H, [l^, it obeys detailed 
balance (shown in Appendix A). Consequently, in the 
T —>■ oo limit, V approaches a unique stationary dis¬ 
tribution, which can be found by applying (N) = 
(N') W(N, N') /IE (N',N) repeatedly. Imposing nor¬ 
malization, we arrive at an explicit, closed form: 

Ni Ne 

7,-(N) = iri(fcd)lJ(p,-!) (9) 

i=i i=i 

where H = E^^jH (fed) LI (pj!) is a ‘partition function.’ 
Note that the particle-hole symmetry (Eq. [7]) is manifest 
here. 

Interpreting P®® as a Boltzmann factor and trivially 
assuming /3 = 1, we can write a ‘Hamiltonian’ 

{ Ni / Ne \ Ne / Ni \ '| 

E In E ! + E In E ! ? 

i=i ) ) 

( 10 ) 

Now, this expression immediately alerts us to the level of 
complexity of this system of ‘Ising spins,’ as H contains 
a peculiar form of long range interactions. Each ‘spin’ 
is coupled to all other ‘spins’ in its row and column, via 
all possible types of ‘multi-spin’ interactions! We are not 
aware of any system in solid state physics with this kind 
of interactions. It is remarkable that such a complex % 
emerges from an extremely simple model of social interac¬ 
tions. Meanwhile, it is understandable that computing H, 
let alone statistical properties of macroscopic quantities, 
will be quite challenging (Appendix B). Nevertheless, as 
the next two sections show, we are able to exploit mean- 
field approaches to gain some insight, and to predict some 
macroscopic observables. For generic {Nj,Ne), we find 
excellent agreement with simulation data. As in standard 
equilibrium statistical systems, these theories fail in the 
neighborhood of critical points, which turn out to be the 
TV/ = Ne line here. 


III. STATISTICS OF A, THE TOTAL NUMBER 
OF CROSS-LINKS 

In the XIE model, the most natural macroscopic 
quantity to study is X, the total number of I-E links. In 
addition, its average {X) can serve as an ‘order parame¬ 
ter,’ since it plays the role of the total magnetisation in 
the Ising model. Though there is no natural temperature¬ 
like variable in our model of social networks, there is a 
natural external-field-like control parameter: h (or A). 
Our main interest in this section is how (A) varies with 
Nj^e- In other words, what is the ‘equation of state’ 
m{h;Af) for the XIE model? If it were like the Ising 
model below criticality, m (Ti = 0±; A? —t oo) = img, 
where 0 < rus < I is the spontaneous magnetisation, ac¬ 
companied by ordinary O ( I/VTE) fluctuations. In the 
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FIG. 2. (Color online) The distributions, P{X), compiled 
with the time traces of X for three cases with {Nj, Ne) near 
and at criticality: (101,99) (green dashed line), (99,101) (blue 
dot-dashed line), and (100,100) (red solid line). 


Thouless effect, a discontinuity would be accompanied 
by anomalously large fluctuations. Here, we will find 
that our m displays an extreme Thouless effect in the 
thermodynamic limit, i.e., m (h = 0±; A/” —>■ oo) = ±1, 
together with extraordinary, 0(1), fluctuations around 
m{h = 0) = 0. 

The preliminary results for N = 200, displayed in Figs. 
1-3 in [l9| gave a hint of these remarkable properties. To 
confirm and to improve on those results, we carry out 
much longer runs: discarding the first 5 x 10^ MCS, tak¬ 
ing measurements every 50 MCS, for up to 10^^ MCS (for 
the Nj^e = 100 case). Using these long traces, X (t), we 
find much more information than just the average (X)] 
we obtain a much more accurate picture for the whole 
steady state distribution P (X): Fig. [51 Not surpris¬ 
ingly, they are sharply peaked and Gaussian-like for the 
off-critical cases (green and blue on line), while the dis¬ 
tribution in the Nj = Ne case (red on line) is essentially 
flat over most of the full range, [0, A/]. The flat plateau in 
P (X) gives the impression of an unbiased random walk 
(bounded by ‘soft’ walls near the extremes of the allowed 
region). In stark contrast, characteristic of co-existence 
in an ordinary first order transition, M [t) for an Ising 
system below criticality spends much of its time hover¬ 
ing around the spontaneous magnetizations, ±Msp, and 
makes rare and short excursions from one to the other. 

Another principal characteristic is metastability or 
hysteresis. Though not shown explicitly here, we ob¬ 
serve neither. When the two nodes ‘change sides,’ i.e., 
(101, 99) —?► (99,101), X/Af simply marches from ~ 15% 
to ~ 85% in ~ 3500 MCS. In other words, on the av¬ 
erage, X changes by about two links per MCS. We also 
considered having (5 = 4 or 6 ‘defectors’ instead of just 
two, in systems with N = 400 and 800. In all cases, 
the average ‘velocity’ is approximately 6 per MCS. Intu¬ 
itively, we may attribute this to the action of the S extra 
E nodes, but it remains to be shown analytically. In all 
respects, there is absolutely no barrier between the two 
extremes of A ! As for criticality itself (100,100), the 



FIG. 3. (Color online) Power spectra, I (w), associated with 
the time traces of X{t). The dashed black line is proportional 
to IjijP. Note that the spectra associated with the two off- 
critical cases (farthest from the dashed line, shown as green 
dots and a blue line) are statistically identical, as expected 
from particle-hole symmetry. Both scales are logio- 


notion of X executing a pure random walk (RW) can be 
further confirmed by studying its power spectrum. With 
T = 2 X 10^ measurements (of runs of 2 x 10® MCS), we 
compute the Fourier transform X {uj) and then average 

over 100 runs to obtain / (w) = ^|A ( 0 ;)]^^ In Fig. [31 we 

show plots of log I vs. logw, as well a straight line (black 
dashed) representing The upper set of data points 
(red on line), associated with Nj = Ne = 100, are statis¬ 
tically consistent with the RW characteristic of The 
cutoff at small oj can be estimated from the finite range 
available to the RW (~ 7000 here). Since AX = ±1 
in each attempt, we can assume the traverse time to be 
about 7000^ = 5 x 10^ attempts, or ~ 2.5 x 10® MCS. 
Given that this value is comparable to 1/10 of our run 
time, it is reasonable to expect deviations from the pure 
uj~'^ as we approach uj ~ 10. By contrast, the power 
spectra of the two off-critical cases (lower set of data, 
green and blue on line) are controlled by some intrinsic 
time scale associated with both the restoration to {X) 
and the fluctuations thereabout. Indeed, this I (w) is en¬ 
tirely consistent with a Lorentzian, i.e., oc 1/ {uP' -I-Wq). 
Given our limited understanding of the dynamics of this 
model, estimating wq is beyond the scope of this work. 

As shown in ref. M , some characteristics of an ex¬ 
treme Thouless effect (presence of a jump, absence of 
metastability and hysteresis, existence of a flat plateau in 
P (X), etc.) can be qualitatively understood in terms of 
a crude mean field approximation. That approach starts 
with the exact 

P(A)^^<5(A,S,,n,,)p- (N). 

{N} 


( 11 ) 
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and replaces every by its average X/Af in the sum, 
resulting in an approximate P (X). Its maximum can be 
used to predict (X) and so, an approximate equation of 
state m {h;Af). Remarkably, at the lowest order in l/AA, 
this approach leads to an extreme Thouless effect [llLfl^ . 
i.e., 

m = sign (h) , (12) 

the absence of metastability and hysteresis when h 
changes sign, as well as a flat plateau in P {X) for h = 0\ 
Keeping the next order in F, we find 

m{h]N) = sign{h) - (13) 

which provides qualitatively good agreement with most 
of the N = 200 data set. Clearly, this mean-field ap¬ 
proach captures some key features of the XIE model, 
even though it is not quantitatively reliable. 

Before we turn to a much better mean field theory, let 
us present the rudiments of a more complete portrait, to¬ 
wards a systematic scaling study. Specifically, we study 
m {h; N) in the transition region, using simulations with 
N up to 3200. Of course, unlike the Ising model, there 
is no natural variable in the social network correspond¬ 
ing to temperature. Nevertheless, the results suggest the 
presence of anomalous power laws and possible data col¬ 
lapse. 

To study the behavior near criticality, we measure the 
average {X) for all possible values of Nj^e which corre¬ 
spond to h < 0.01, with N = 200,400, ...,3200. In other 
words, we use the appropriate values of A = 2,4, ...32. 
Starting with a half filled set of I-E links, each point was 
calculated with a run of up to 5 x 10® MCS with mea¬ 
surements made about every iV/4 MCS. Verifying that 
our data is consistent with particle-hole symmetry, we 
present results only for h > 0 and m = 1. Fig. 0^1 
shows results within the small regime h G [0,0.01] and 
m S [0.7,1.0]. We see that, indeed, m approaches 

-1-1 a.s N ^ oo. Note that, we can access smaller h 
in systems with larger iV, as its minimum is 2/N . 

With this set of data, we make two log-log plots (Fig. 
0 }d) showing that (i) for fixed h = 0.01, m —>■ 1 with 
jy-o.yi Jqj. fixed A = 2, m —>■ 1 with The 

solid lines represent linear fits with correlation coefficient 
(R®) greater than 0.999. Using this information, we plot 
against in Fig. 0}:. Though somewhat 

rough, this plot does qualitatively indicate data collapse 
and provides the first steps towards an in-depth finite 
size scaling analysis. Such a study is beyond the scope 
of this paper and will be reported elsewhere [l^. If the 
exponents found here are confirmed, they signal a signif¬ 
icant deviation from the mean-field values (e.g., Ea fT3l) : 
(0.70,0.36) as opposed to (1,0). Of course, they also 
provide fertile ground for renormalization group analy¬ 
sis, along the lines of ref. 0[il. 



FIG. 4. (Color online) (a) m{h) for various N’s and allowed 
A’s with h < 0.01. (b) Log-log plots of indicated subset of 
points in (a). Solid and dashed lines are linear fits with slope 
—0.363 and —0.707, respectively. Both correlation coefficients 
are greater than 0.999. (c) Scaled plots of the points in (a), 
showing tolerable data collapse. Dashed line is linear, as guide 
to the eye. 


IV. DEGREE DISTRIBUTIONS: SIMULATION 
RESULTS AND A DYNAMIC MEAN-FIELD 
THEORY 

In this section, we turn to a ‘mesoscopic’ quantity 
which offers a more detailed perspective than the macro¬ 
scopic X, as well as major contrasts between Ising-like 
statistical systems and those associated with networks. 
For the Ising model, a natural quantity to study is the 
total magnetisation, which corresponds to X. But, it 
is not usual to study the statistics of the magnetisation 
in a row or a column. Yet, for the XIE model, the 
corresponding quantity is the degree distribution, p(fc), 
which is one of the most common ways to characterise 
a network. Thus, we devote the rest of this paper to 
these distributions, illustrating with simulation data (for 
N = 200), as well as offering a more effective mean-field 
theory. Specifically, unlike the mean-field approach for 
P {X), we will formulate an approximate dynamics for p 
and arrive at much better agreement with data, for all 
Nj ^ Ne- In particular, m (0.01; 200) above differs from 
data by ~ 15%, while the theory below produces a value 
within 0.02%. 

Before presenting the data, let us set the stage for dis¬ 
cussing two degree distributions. In general, associated 
with a network with several subgroups or communities, 
we can study many such distributions, to describe the 
various intra- and inter-community links. For X/U, the 
intra-community links are static and so, we need to study 
only two: p/(fc/) and pE{kE), related to the degrees of the 
/’s and U’s, respectively. To illustrate, a network with 
pi{l) = 5 and pe{^) = 1 is shown in Fig.[TJ From the p’s, 
average degrees (fc/) and (/c^;) can be found. Note that 
{ki) ^ (fcs) typically, but they are related, in the steady 
state, by the following. Since TV/ (fc/) is just the average 
number of cross-links, while an extrovert has fc/; — Ne +1 
links to the I’s, we have 

TV/ (fc/) = {X) = Ne (kE) — Ne {Ne — 1) (14) 
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This complication can be bypassed, especially in view of 
the particle-hole symmetry discussed above, by introduc¬ 
ing a ‘hole’ distribution for the E^s: Ce (pe)- Clearly, ( 
is intimately related to p, namely, C{p) = P ~ ^ ~ p)- 
Meanwhile, Ne (pe) = M — (X), so that a manifestly 
symmetric constraint on p, is Nj (kj) + Ne (pe) = A/”. 

Next, let us present the simulation results. Starting 
with a network with various initial conditions (null graph, 
complete graph, random half-filled), we evolve the system 
according to the simple rules given above. Not surpris¬ 
ingly, after O {N) MCS, all the /-/ links are absent while 
all the E-E links are present. To be quite certain that the 
system has equilibrated, we discard the first 5 x lO”^ MCS. 
Thereafter, we measure the degrees of each node every 50 
MCS. The distributions are then obtained as the average 
over 10® measurements. Shown in Fig. EKa) are p (k) for 
three cases: {Ni,Ne) = (150,50), (125,75), (101,99). 
Evidently, each p consists of two components, associated 
with k < Ne and k > Ne — 1- Generally, these com¬ 
ponent are disjoint and so, they can be identified unam¬ 
biguously with respectively, pi and pE (shown with open 
and solid symbols). Note that, in the (101,99) case, de¬ 
spite having just two members less, the ~ 100 extroverts 
are unable to create enough cross-links, so that there are 
essentially no I’s with say, 50 links! Of course, this re¬ 
sult is entirely consistent with our observations above, 
showing that a typical I has only about 15 links. Apart 
from having these two components, the most prominent 
feature is that neither component resembles the 

degree distribution of a homogeneous population with 
preferred degree k [mil- As will be shown, they are 
well approximated by Poisson distributions, an analytic 
result of our mean-field theory. 

What happens when the introverts ‘defect’? Though 
the changes appear dramatic, they should not be a sur¬ 
prise, given the underlying particle-hole symmetry in 
XIE. Thus, we illustrate in Fig. EKb) the degree dis¬ 
tribution for {Ni,Ne) = (99,101) (blue circles), as well 
as the previous case of (101,99) (green circles). By ex¬ 
changing Nj ^ Ne and plotting the degree distribution 
vs. p = N — 1 — k, we find perfect (within statistical 
errors) overlap between the blue and green data points. 


A. Self-consistent mean-field approximation 
(SCMF) 

Given the exact steady state distribution (Eq. E]), the 
p’s can be computed, in principle, from e.g., pi (kj) = 
^ ~ (^) *)• practice, this 

task is as difficult as computing P (X), so that we again 
resort to a mean-field approach. The main difference be¬ 
tween the earlier scheme and this one is that the approxi¬ 
mation will be applied to the underlying dynamics of the 
model [H, [l^ (as opposed to evaluating the sum 
above). In other words, we formulate an approximation 
on the transition probabilities - for the degree of a partic¬ 
ular node to increase/decrease by unity: R{k —>■ fc± 1). 


Once these are determined, we impose the steady state 
condition 

p{k) R{k ^ k — 1) = p{k — 1) R{k — 1 ^ k) (15) 

to find p (fc) (being an approximate p, again denoted by a 
tilde) in closed form. The strategy is as follows. Exploit¬ 
ing particle-hole symmetry, we will consider the two dis¬ 
tributions, Pi and C,e, as well as two sets of rates, Ri^e- 
Each R will depend on an unknown parameter, repre¬ 
senting the average degree of the opposite community. 
From these, explicit expressions for pi and (^e can be ob¬ 
tained. In turn, the average degrees can be computed 
and the unknown parameters can now be fixed through 
self-consistency. In this spirit, we refer to this scheme as 
a SCMF approximation, details of which can be found in 
Appendix C. Here, we simply quote the results: 

^Ne—^j ^Nj—PE 

Z, (Ns:-Zs(N,-psy. 

(16) 

where A, p are constants which can be obtained from 
Nj^e alone and Z’s are normalization factors (Eqs. 
IC4IC10I) . Both are truncated Poisson distributions, since 
ki £ [0, Ne] and pe £ [0, ?V/] • Instead of quoting A and p 
from the SCMF calculation, we plot the full distributions 
predicted by Eqs. (IC5IC9I) . shown as solid black lines in 
Fig. EKa). We should emphasize that no fit parameters 
have been introduced in this approach; the lines depend 
only on the control parameters, Nj^e- It is clear that the 
agreement between theory and simulation data is excel¬ 
lent for Nj > Ne- By symmetry, it will also be quite 
good for cases with Nj < Ne- Indeed, disagreement be¬ 
tween theory and data is visibly detectable only in the 
tails of the next-to-critical case, (101,99), a sign that cor¬ 
relations can no longer be entirely neglected here. From 
these distributions, we easily obtain {X) using Eq. (1141) . 
Unlike the results from the previous Section, all of these 
predictions fall within statistical errors of the data. Since 
mean field schemes are not the start of a systematic set 
of approximations, it is unclear why the approach here 
is so much more successful at capturing the essentials of 
the model. 

To end this section, let us illustrate, with the sym¬ 
metric case (100,100), the challenges of ‘criticality.’ As 
expected, there are non-trivial obstacles for both Monte 
Garlo simulations and theoretical understanding. First 
of all, it takes much longer for the system to settle, typ¬ 
ically a hundred-fold longer than the Nj Ne cases. 
To compile a reliable histogram for the p{k), shown in 
Fig. El we take 10^° measurements in a combination of 5 
runs, each of which lasts for 10^^ MGS (after discarding 
the first 10^ MGS for the system to settle into steady 
states). Note that, to untangle pi^E in the central pair 
of points, we recorded separately whether an / or a A 
node has 99 or 100 links. Of course, the distribution is 
symmetric, i.e., pi = Ce- Unlike the off-critical cases, 
these distributions display broad and flat plateaux. Un¬ 
doubtedly, the physics underlying these also gives rise 










FIG. 5. (Color online) Degree distributions, p, for several cases with Nj + Ne = 200. Simulation results for the low/high 
k components, associated with introverts/extroverts, are denoted by open/solid symbols, (a) The symbols for {Ni,Ne) are 
orange triangles (150,50), purple squares (125,75), and green circles (101,99). The solid black lines are predictions from a 
self-consistent mean-field theory, (b) When two introverts ‘change sides,’ a dramatic jump in p{k) results, with the case of 
(99,101) shown as blue diamonds. 



FIG. 6. (Color online) Simulation results of the degree dis¬ 
tribution for the symmetric (100, 100) case. The two compo¬ 
nents, denoted by open and solid diamonds, can be associated 
with the separate distributions pi and pE, respectively. 


to the plateau in P{X). Quantitative understanding of 
these large fluctuations remains elusive, while the SCMF 
prediction for this case is, not surprisingly, far from ideal. 


V. SUMMARY AND OUTLOOK 

In this article, we report findings concerning an ex¬ 
traordinary phase transition in a minimal model of dy¬ 
namic social networks with preferred degree. Consisting 
of extreme introverts/extroverts, who only cut/add con¬ 
nections, these dynamic networks have only the I-E links 
and reduce to bipartite graphs. With only two control 
parameters, this seemingly trivial model displays 

surprising behavior. In particular, we find compelling 
evidence that, in the limit of large populations, (i) the 


likelihood of a link being present, (x), jumps discontin- 
uously, from 0 to I, when Nj drops below Ne and (ii) 
in the Nj = Ne case, x assumes all values in [0,1] with 
equal probability. Such remarkable properties have been 
observed in other statistical systems, e.g^MD Ising mod¬ 
els with certain long range interactions . We can 

place the similarity between our dynamic XIE model 
and equilibrium Ising systems on firmer grounds, by map¬ 
ping their microscopic configurations one-to-one and re¬ 
garding the evolution of the former as Glauber spin-flip 
dynamics on the latter. Thanks to restoration of detailed 
balance, we are able to find an exact expression for the 
XIE stationary distribution, Interpreting it as a 

Boltzmann factor, — InP®'* can be regarded as a ‘Hamil¬ 
tonian’ for the spins, though much more complex than 
typical Ising models. Statistical properties of our system 
can now be explored along standard routes. 

To study of the collective behavior arising from such 
microscopies, we focus on the degree distributions and x, 
the fraction of cross-connections. Since x corresponds to 
m (the magnetisation in the Ising model), while hor A 
(Eq. 0]) can serve as an external magnetic field, a nat¬ 
ural question for us is: What is the ‘equation of state,’ 
m{h)? Though naive expectations lead to the trivial 
m = h, both simulations and mean field theories point 
towards the contrary: m = sign ( /i), which is a hall¬ 
mark of an extreme Thouless effect [ll|, [l^ • Apart from 
the macroscopic X, we also study degree distributions 
p(k), ‘mesoscopic’ quantities which are commonly used 
in characterizing networks. Remarkably, the predictions 
from a mean field approximation, formulated at the level 
of the underlying dynamics for p(fc,t), are in excellent 
agreement with data (for all Nj ^ Ne)- 

The results of this first study are encouraging and pro¬ 
vide us with good stepping stones towards more system¬ 
atic investigations. The most obvious question may be 





what forms do thermodynamic limits take, assuming they 
exist. Preliminary studies [2^ suggest that, for fixed h, 
the degree distributions (11611 approach non-trivial lim¬ 
its. On the other hand, if A were held fixed instead, it 
is unclear what the limiting behavior is. How the sys¬ 
tem approach such limits is the next issue. A detailed 
study of finite size scaling should be undertaken , us¬ 
ing both simulations and theoretical techniques. Pursu¬ 
ing an exact computation of the partition function, and 
perhaps P{X), poses a worthy challenge. The failure 
of mean-field theory, especially near ‘criticality’ hints at 
the importance of correlations. However, there is no spa¬ 
tial structure in our model and so, the usual notion of 
correlation length is ill defined. Nevertheless, we have 
some evidence of strong correlations, in the sense that 
the joint distribution, p{kj,kE), can be quite different 
from the product pi (kj) pE (kE)- Systematic investiga¬ 
tions of them are straightforward and worthwhile. From 
a theoretical point of view, it would be desirable to de¬ 
velop an understanding for why the SCMF is so much 
more successful than the standard mean-field approach 
(Section HI). Such insights may have impact beyond this 
study, as they may reveal how best to formulate mean- 
field approximations. 

Beyond exploring these questions, we can extend the 
XIE model in an orthogonal direction, arguably of 
purely theoretical interest (at present). We may treat 
TL as a genuine Hamiltonian in a standard study of crit¬ 
ical phenomena in thermal equilibrium. In other words, 
we propose to study the statistical mechanics of a L x L 
system associated with the Boltzmann factor 

V ocexp{-P[n- BX]} (17) 

Here, P is the usual inverse temperature variable, while 
the bias B plays the role of a symmetry breaking, ‘mag¬ 
netic field’ (similar, but not identical to h in the XIE). 
It is interesting to note that, while the critical control pa¬ 
rameters of a typical system (e.g., Tc in Ising, Tc and Pc 
for liquid gas) are not known, they are given precisely by 
/3c = 1 and B^ = 0 here. For this ‘purely theoretical’ sys¬ 
tem, work is in progress , to explore the usual avenues 
of interest: static and dynamic critical exponents, scaling 
functions, universality and the classes, etc. In the lan¬ 
guage of renormalization group analyses (which proved 
to be highly effective in dealing with other mixed-order 
transitions ESi), we already know that B lies on the 
critical sheet and can inquire about fixed points and their 
neighborhoods, irrelevant and relevant variables (e.g., if 
there are others besides — 1 and B), etc. In addition, 
there are unusual challenges, such as the lack of a natural 
correlation length in such a system. 

Beyond the XIE model and its purely theoretical com¬ 
panion, there is a wide vista involving dynamic networks 
with preferred degrees. For instance, instead of assign¬ 
ing one or two k’s to a population [l^, , it is more 

natural to assign a distribution of k’s. There are also 
multiple ways to model interactions between the vari¬ 
ous groups. For example, even with just two groups, it 


is realistic to believe that an individual may have two 
preferred degrees, one for contacts within the group and 
another for those outside. Surely, this kind of differen¬ 
tial preference underlies the formation of social cliques. 
Beyond understanding the topology and dynamics of in¬ 
teracting networks of the types described here, the next 
natural step is to take into account the freedom associ¬ 
ated with the nodes, e.g., opinion, wealth, health, etc., on 
the way to the ambitious goal of understanding adaptive, 
co-evolving, interdependent networks in general. Along 
the way, we can expect the unexpected, such as the emer¬ 
gence of the extreme Thouless effect in this XIE model, 
arguably the simplest of all interacting social networks. 


Appendix A: Restoration of detailed balance 

In this appendix, we show that all Kolmogorov loops 
are reversible in the XIE model and so, detailed bal¬ 
ance is restored [T7I| . Since the full dynamics occurs 
on the Af cross-links, the configuration space consists 
of the corners of an A/'-dimensional unit cube, while 
adding/cutting a link is associated with traversing an 
edge therein. Clearly, products of the ratios of forward 
and reversed transition rates around any closed loop can 
be expressed in terms of those around ‘elementary loops’ 
- i.e., loops around a plaquette on the A/’-cube. We will 
show that the ratio associated with every plaquette is 
unity and so, all Kolmogorov loops are reversible. 

First, it is easy to see that if an elementary loop con¬ 
sists of modifying two links connected to four different 
nodes, then the actions on each link are unaffected by 
the other. In other words, rates associated with oppo¬ 
site sides of the square (loop) are the same. Thus, their 
product in one direction is necessarily the same as in the 
reverse. We need to focus only on situations where the 
two links are connected to three nodes, e.g., ij and im. 
For any such loop, let us start with a configuration in 
which both are absent (riy = riim = 0). Let the states 
of node be such that i has ki links, and j and m have pj 
and Pm ‘holes’, respectively. Then one way around the 
loop is adding these two links followed by cutting them, 
which can be denoted as the sequence 



and leaving the rest of N unchanged. The associated 
product of the transition rates is, apart from an overall 
factor of 


Pj Pm ki + 2 ki + 1 

Now, the reversed loop can be denoted as 




associated with the product 

111 1 
Pm Pj ki + 2 ki + 1 


(A4) 


is given by Q (z) /H, where 
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which is exactly equal to Eqn. (IA2I) . From symmetry, 
we can expect the same results for loops involving two 
introverts and one extrovert (i.e., ij and kj). Thus, we 
conclude that the Kolmogorov criterion is satisfied and 
detailed balance is restored in this XIE limit. Our system 
should settle into a stationary distribution without prob¬ 
ability currents, much like the Boltzmann distribution for 
a system in thermal equilibrium. 


Appendix B: Considerations for computing O and 

P{X) 


Ni Ne 

(2;) = ^ n (^*0 n 

{ny} i=l i=l 


I.e., 


VuDv exp < — 


id 


1 / N 


(B5) 

Of course, this integral is fraught with the same issues as 
in n. Obviously, even if G (z) can be found, inverting it 
may pose other challenges. Thus, deriving the presence 
of a plateau in P{X), as well as anomalous exponents 
associated with its edges, will be a non-trivial endeavour. 


Exploiting 


k\= due-'^u*^ 

Jo 


(Bl) 


the ‘partition function’ can be expressed as 


ft — / dUiC 


n 


dvjC 


{nij} i 


where we have used Eqn. (|6l). Exchanging the configu¬ 
ration sum and integrals, we can perform the former to 
find 


n = 


/ 


VuVve 


-{T,iUi + T.j Vj] 


En 

{ny} 


which can be cast as 


VuVv exp < — 


E 


^7 1 / \ 


Here, 


'DuDv = / dui n dvj 


(B2) 


is a precursor for functional integrals if we take the con¬ 
tinuum limit Ui ,Vj -d u (x ), v (y) and regard the resul¬ 
tant as a two component, 1-D field theory. An attempt to 
use standard steepest descent leads to the following com¬ 
plications. As long as Nj 7 ^ Ne, the maximum of the 
integrand in (|B|) is located at a boundary of the region of 
integration. However, for Nj = Ne, the maximum is a 
line given by Ui = u, Vj = v and u + v = Nj^e- Both are 
non-standard behavior and require more care to proceed. 

Similar considerations can be given to the computation 
of P (A), in the sense that its generating function 


G{z) = Y,z^P{X) (B3) 


Appendix C: Degree distributions in a self 
consistent mean-field approach 


In this appendix, we provide some technical details for 
the SCMF scheme. Note that parameters (A, p) in the 
main text are, in fact, the physically meaningful quanti¬ 
ties {{peY , {kiY) here. 

First, consider a particular I node, with pi (ki) being 
the probability to find it having kj links. Then, provided 
ki > 0, Ri {kj ^ ki — 1) = 1/N which is the probability 
that this node is chosen to act. By contrast, the exact 
rate for having a link added (fc/ — 1 —>• kj) is more com¬ 
plicated, since it depends not only on all the Ne — kj + 1 
extroverts not connected to it, but also on how many 
‘holes’ each has - through l/pj (in Eq. [5]). To proceed, 
we make judicious approximations. In the spirit of mean- 
field theory, we can replace 1/pj by the average {I/peY, 
where the prime stands for an average restricted to nodes 
with pe > 0. Though we can formulate the theory with 
{1/pe)', let us make a further simplifying approximation 
and replace it by 1/ (p£;)^ So, we write 

1 (Cl) 

A (pE) 


So far, {peY is an unknown parameter. If we had the 
distribution of an extrovert’s holes, (^e [pe), then we have 
the following relation: 


/ y _T.pe>oPe(e{pe) 
" J:pe>oCe{pe) 


(pe) 

1 - Cs ( 0 ) 


(C2) 


But, Ce (pe) is unknown. Nevertheless, at this stage, we 
can exploit Eq. (1151) and readily find 


PI {ki) 


Ne — kj + 1 Ne — kj + 2 

// ,r\NE — ki 

[{pe) ) 

[NE-ki)l 


Ne 


PI (0) 


oc 


(C3) 
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Before continuing to study Re, let us work with this 
expression further. Since pi = Ne — fc/ is the number 
of ‘holes’ associated with an I node, we recognize this 
as a Poisson distribution (truncated aX Ne) for the hole 
distribution. Imposing normalization, we find a compact 
closed form, Q {pi) = {{peYY^ /Zipil, where 

zi = j2i-,{{py)'/i'- (C4) 

is the sum of the first Ne + 1 terms of an exponential 
series. Despite its simplicity, the notation Q {pi) may be 
too confusing and so, we will quote the final result for pi 
as 


Pi {ki) 


{, ,/\NE—ki 

[{pe) ) 

Zi{NE-ki)\ 


(C5) 


with {pe)' being a to-be-determined parameter. 

Next, we turn to a particular E node and, exploiting 
‘particle-hole’ symmetry, consider its hole distribution, 
Cb {pe)- Since adding a link is decreasing pE by unity, 
we again have Re {pe + 1 — t Pe) = ^/N, the probability 
that this node is chosen to act, provided pE > 0- Mean¬ 
while, it is connected to Ni — pE (i-e., kE — Ne + 1) 
introverts, each of which has ki links. As above, we rely 
on the same arguments and replace the ki’s by a suitable 
average: 


Re {pe — t Pb + 1) 


Ni -pE 1 

N {kiY 


where 


{kiY 


(ki) 

1 - PI (0) 


Recasting Eq. m for we have 

Cb (pb) = TT——Cb (pb + 1) 

Ni - PE 


(C6) 


(C7) 


(C8) 


Again, this recursion relation leads to a (truncated) Pois¬ 
son distribution in Ni —pE, and imposing normalization, 
we have explicitly 


Cb (pb) 


{{kiYf’~^^ 

Ze {Ni — peY- 


(C9) 


with 

ZE = J2flo{ikYY/i'. (CIO) 


Of course, (ki)' here is also an unknown, to-be- 
determined, parameter. Note that, along with Eq. (IC5() . 
this result again confirms the underlying particle-hole 
symmetry. 

Finally, we make the last approximation. Instead of 
the exact (and unknown) parameters, (jje)' and (ki)', let 
us approximate them by using pi and Cb in Eqs. ()C7IC2ll 
instead. Since p/ and Cb depend on {peY and (ki)', 
respectively, we may define the functions / and g: 


{kiY 

{peY 


Ekipi {ki) _ 

I-pz(O) = 
Z>peCe {pe) 
1 - Cb (0) 


/ ((pb)') 

= 9i{kiY) 


(Cll) 

(C12) 


Making a plot of these functions in the {k)'-{pEY 
plane, the point of intersection then determines, self- 
consistently, the values for these two parameters. In 
practice, it is simple to start with, say, a trial value po for 
{peY and compute (ki)' through Eq. (IC5D . Inserting this 
(ki) into Eqn. (IC9I) . we compute Cb and the associated 
{peY■ If this result is not po, then vary the latter until 
they agree. In other words, this process will lead us to 
the solution: {peY = p(/((pb)'))- Substituting these 
values {{peY and (ki)') into Eqs. (IC5IC9I) . the degree 
distributions can be plotted. 
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